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ABSTRACT 

We analyze 24 type I X-ray bursts from GS 1826—24 observed by the Rossi X-ray Timing Explorer 
between 1997 November and 2002 July. The bursts observed between 1997-98 were consistent with 
a stable recurrence time of 5.74 ± 0.13 hr. The persistent intensity of GS 1826—24 increased by 36% 
between 1997-2000, by which time the burst interval had decreased to 4.10 ± 0.08 hr. In 2002 July the 
recurrence time was shorter again, at 3.56 ± 0.03 hr. The bursts within each epoch had remarkably 
identical lightcurves over the full 150 s burst duration; both the initial decay timescale from the 
peak, and the burst fiuence, increased slightly with the rise in persistent flux. The decrease in the burst 
recurrence time was proportional to M ~ 1 - 05±0 - 02 (assuming that M is linearly proportional to the X-ray 
flux), so that the ratio a between the integrated persistent and burst fluxes was inversely correlated with 
M . The average value of a was 41.7 ± 1.6. Both the a value, and the long burst durations indicate that 
CN \ the hydrogen is burning during the burst via the rapid-proton (rp) process. The variation in a with M 

implies that hydrogen is burning stably between bursts, requiring solar metallicity (Z ~ 0.02) in the 
accreted layer. We show that solar metallicity ignition models naturally reproduce the observed burst 
energies, but do not match the observed variations in recurrence time and burst fiuence. Low metallicity 
models [Z ~ 0.001) reproduce the observed trends in recurrence time and fiuence, but are ruled out by 
the variation in a. We discuss possible explanations, including extra heating between bursts, or that the 
OO ' fraction of the neutron star covered by the accreted fuel increases with M . 

Subject headings: stars: individual (Ginga 1826—238, GS 1826—24) — X-rays: bursts 
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1997), although this source may also have been the origin 
of X-ray bursts observed much earlier by OSO-8 (Becker 
et al. 1976). Optical photometry of the V = 19 coun- 
terpart revealed a 2.1 hr modulation, as well as optical 
bursts (Homer et al. 1998). The delay time measured be- 
tween the X-ray and optical bursts is consistent with the 
binary separation for a 2.1 hr orbit (see also Kong et al. 
2000). Based on optical measurements, the distance, d, 
to the source is at least 4 kpc (Barret et al. 1995). An 
upper limit of 8 kpc has been derived from the peak fluxes 
of bursts measured by BeppoSAX, ASCA and RXTE (in 
't Zand et al. 1999; Kong et al. 2000), which would po- 
sition the source just outside the Galactic bulge. The 
global accretion rate inferred from the observed flux is then 
M w lO^Moyr-^Fp/Sx 10~ 9 ergs cm" 2 s" 1 )^ kpc) 2 , 
about a factor of ten higher than the expected time- 
averaged accretion rate for a NS orbiting a main sequence 
star in a 2.1 hr orbit (Verbunt & van den Heuvel 1995). 
Such a discrepancy is not surprising for a transient accre- 
tor. 

Extensive BeppoSAX Wide Field Camera (WFC) obser- 
vations between 1996 August and 1998 October revealed 
around 70 bursts, which recurred on an average interval 
of 5.76 ± 0.62 hr (leading to the pseudonym "the Clocked 
burster" ; Ubertini et al. 1999). Analysis of the entire WFC 
sample of 260 bursts revealed that the source consistently 
exhibits approximately periodic bursts with a recurrence 
time which decreases significantly as the persistent flux in- 
creases (Cornelisse et al. 2003). Such regularity of bursts 
over a long time span is not observed in other sources 
and implies a remarkably stable accretion rate, as well as 
complete burning of the accreted fuel and a constant fuel 
covering fraction. 

Based on the decay time scale of the bursts observed 
by RXTE (Kong et al. 2000) and the recurrence times, 
Bildsten (2000) suggested that the bursts must involve 
mixed hydrogen/helium (H/He) burning. In this scenario, 
the initial helium flash generates seed nuclei for hydro- 
gen burning, which gives rise to an additional energy re- 
lease and also allows an active rapid proton (rp) process 
of Wallace & Woosley (1981). This process burns hydro- 
gen by successive proton capture reactions on seed nuclei 
that subsequently /3-decay. These reactions move the seed 
nuclei up the proton-rich side of the periodic table (much 
like the r-process which occurs by neutron captures on the 
neutron-rich side of the periodic table). Theoretical work 
indicates that the end point of the rp-process is set either 
by the complete consumption of the available hydrogen 
(Hanawa et al. 1983; Rembges et al. 1997; Schatz et al. 
1998; Koike et al. 1999) or by reaching the closed SnSbTe 
cycle found by Schatz et al. (2001). Uncertainties still re- 
main in the actual total time to burn the matter, since 
the /3-decay lifetimes of many of these proton-rich nuclei 
are not known experimentally. Thus, GS 1826—24 is an 
ideal candidate to test, and provide input for, our current 
theoretical understanding of the role of the rp-process in 
thermonuclear X-ray bursts. 

We present new measurements of the burst recurrence 
time in GS 1826—24 with the Rossi X-ray Timing Ex- 
plorer (RXTE). We first compare the properties of differ- 
ent bursts to highlight their unusually uniform properties. 
We examine the long-term trend of burst recurrence time 



as a function of persistent flux, and test for variations in 
the a-value within different epochs. By comparing to the- 
oretical models of type I bursts, we are able to constrain 
the metallicity of the accreting material, and the fraction 
of the neutron star's surface over which accretion takes 
place. 



2. OBSERVATIONS 

We analysed observations of GS 1826—24 by the All 
Sky Monitor (ASM; Levine et al. 1996), the Proportional 
Counter Array (PCA; Jahoda et al. 1996) and the High- 
energy X-ray Timing Experiment (HEXTE; Rothschild 
et al. 1998) aboard RXTE. The ASM consists of three 
Scanning Shadow Cameras (SSCs) sensitive to 2-10 keV 
photons mounted on a rotating platform, which makes 90 s 
observations ( "dwells" ) covering most of the sky every few 
hours. Data from each SSC from each dwell are averaged 
to obtain the daily intensities of all known sources in the 
field of view. The PCA is made up of 5 identical co-aligned 
Proportional Counter Units (PCUs) sensitive to photons 
in the energy range 2-60 keV, and with a total area of 
i=a 6000 cm 2 . Data is collected during pointed observations 
in standard data modes (1 and 2), as well as a variety of 
user defined modes offering time resolution up to 1 fis and 
256 energy channels. The HEXTE comprises two clusters 
with 1600 cm 2 total area, each of which contains 4 scin- 
tillation detectors collimated to view a common direction 
and sensitive to photons in the range 15-250 keV. 

We searched for X-ray bursts in 1-s binned "Standard- 
1" mode data from 43 public RXTE/PCA observations of 
GS 1826-24 between 1997 November 5 and 2000 Septem- 
ber 27, totalling 347 ks. We also searched in a 30 ks pro- 
prietary observation from 2002 July 27 (observation ID 
70044, PI: Rothschild). We identified 24 bursts in total. 
In three cases the PCUs were switched on during what 
appeared to be the decay from a burst. We estimated 
the start time of these burst candidates by comparing the 
decay light curve with the profile from the burst imme- 
diately before or after. While no blackbody temperature 
decrease was observed in any of these three cases, no de- 
crease was observed for the other bursts when at the same 
stage. Since the observed profiles of the other bursts are 
so uniform, we expect the inferred start times to be cor- 
rect to within 5-10 s. In a fourth case the burst began 
while the satellite was still slewing to the source. For that 
burst the measurement of the peak flux and fluence are 
suspect, but the start time was precisely measured. For 
the analysis presented here, we include these four bursts 
for calculations of recurrence times but exclude them from 
any analysis involving measurements of the flux or fluence. 

For each burst where we observed the entire profile we 
extracted PCA spectra binned initially at 0.25 s (with 
progressively larger bins towards the tail of the burst) 
from Event mode data (125 [is time resolution and 64 en- 
ergy channels). We calculated a separate response ma- 
trix for each burst using PCARSP version 8.0 (supplied 
with lheasoft version 5.2, 2002 June 25). We fitted 
an absorbed blackbody model to each spectrum, using 
a persistent spectrum extracted from a 16-s interval be- 
fore the burst as background. This approach is relatively 
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Fig. 1.— Profiles of 20 X-ray bursts from GS 1826-24 observed 
by RXTE between 1997-2002, plotted with varying vertical offsets 
for clarity. The upper group of 7 bursts were observed in 1997-98, 
the middle group of 10 bursts in 2000, while the lower group of 3 
were observed in 2002. The bursts from each epoch have been time- 
aligned by cross-correlating the first 8 seconds of the burst. Error 
bars indicate the lcr uncertainties. 

standard for X-ray burst analysis (e.g. Kuulk- 
ers et al. 2002). We estimated the bolometric 
burst flux at each timestep as Fb u = 1-0763 x 
10~ n T^ hi Kbb,i crgscm~ 2 s _1 , from the blackbody 
(color) temperature Tbb,i and normalization Kbbd = 
(i?bb/rfiok pc ) 2 ) where Rbb is the apparent photosphere 
radius in km, and (iiokpc the distance to the source in 
units of 10 kpc. We define the burst start as the time at 
which the bolometric flux exceeds 25% of the peak flux, 
and the rise time as the interval for the flux to subse- 
quently exceed 90% of the maximum. We fitted the flux 
decay to an exponential profile with a "break" at which 
point the e-folding time is allowed to vary discontinuously. 
We integrated the measured fluxes over 150 s cover- 
ing the burst, and extrapolated beyond this based on the 
exponential fits, to derive the fluence. 

We also estimated the instrumental background us- 
ing PCABACKEST version 3.0 and the "combined" bright 
source models, and measured the (absorbed) persistent 
2.5-25 keV PCA flux by integrating over an absorbed 
blackbody plus power law model fitted to the persistent 
(pre-burst) spectra. The mean reduced-% 2 for the persis- 
tent spectral fits was 1.07 (56 degrees of freedom). The 
neutral column density was in most cases poorly con- 



strained and not significantly different from zero, and in 
the mean was nn = (2.4 ± 1.4) x 10 22 cm~ 2 . While this 
model provided a good fit to the PCA data alone, com- 
bined fits including the HEXTE spectrum additionally re- 
quired modelling of the high-energy spectral cutoff (see 
section §3.2). 

3. BURST PROFILES, ENERGETICS AND RECURRENCE 

TIMES 

The X-ray bursts observed by RXTE were remarkably 
similar to each other (Fig. 1). The rise times were rel- 
atively long, betseen 4.75 and 7 s (5.6 ± 0.6 s on av- 
erage). The first exponential decay timescale increased 
from 14.7 ± 0.7 to 17.5 ± 1.1 s between the 1997-98 and 
2000 bursts, and to 19.1 ± 1.3 for the 2002 bursts. The 
variation of the burst profile with epoch is obvious in 
the averaged lightcurves (Fig. 2). The second expo- 
nential timescale was, on average, 43 ± 1 s. The peak 
fluxes also showed weak evidence for a decrease with 
time; the mean for the 7 bursts observed in 1997-98 was 
(33.0 ± 0.8) x 10~ 9 ergscm~ 2 s _1 , while for the bursts ob- 
served in 2000-2 it was (30.5 ± 1.1) x 10~ 9 crgscm^s- 1 
(note that the averages of burst properties calculated here 
exclude the bursts which we did not observe in their en- 
tirety). This decrease was substantially larger than the 
variation in the pre-burst persistent emission (see §3.1, be- 
low). Thus, it appears unlikely that the observed variation 
in the peak burst flux arose as a side-effect of subtracting 
the persistent emission as background. The net effect of 
the variations in peak flux and timescale was to keep the 
fluence approximately constant, at ~ 1.1 x 10~ 6 ergs cm~ 2 . 
None of the bursts exhibited evidence for radius expan- 
sion, so that the maximum burst flux is a lower limit to 
the Eddington luminosity. The implied distance limit is 
consistent with that derived from previous observations. 
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Fig. 2.— Mean profiles of 7 X-ray bursts from GS 1826-24 ob- 
served by RXTE during 1997-98 (grey histogram), and of 10 bursts 
observed during 2000 (black histogram). The bursts from 2002 have 
similar profiles to those from 2000. Error bars indicate the lcr un- 
certainties, derived from the scatter of the flux within each time bin 
over all the bursts. The inset shows the same profiles, expanded to 
show more detail around the burst rise and peak. 
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Fig. 3. — (top panel) Recurrence time for thermonuclear bursts 
from GS 1826-24. Shown are the 1997-2002 measurements by 
RXTE (open squares), as well as earlier measurements by Bep- 
poSAX (in 't Zand et al. 1999; Cocchi et al. 2001) and ASCA (Kong 
et al. 2000). (bottom panel) Long-term ASM count rate history of 
GS 1826—24. 10-day averages were calculated from 1-day measure- 
ments, excluding those points with errors > 0.7 counts . The 
dot-dashed line shows the linear fit curve. The linear trend is sig- 
nificant at the lOcr level. 



3.1. Long-term burst interval history 

Due to its low Earth orbit, RXTE can typically only ob- 
serve any given source for 65% of each 90 min orbit. Thus, 
it is likely that some bursts were missed during the gaps 
between observations. The shortest burst intervals were 
found between pairs of bursts observed on 1997 Novem- 
ber 5-6, 2000 Jun 30, and 2002 July 29, at At = 5.88, 
4.00 hr, and 3.58 hr, respectively. While the first value 
is consistent with other measurements around the same 
time (e.g. Ubertini et al. 1999; Cocchi et al. 2001), the 
latter two are almost a third shorter. The other burst in- 
tervals measured from the RXTE observations are at least 
a factor of 2 greater than the shortest intervals, and were 
close to integer multiples of the two shortest intervals in 
each epoch. Furthermore, in each of the longer burst in- 
tervals, the predicted intermediate burst times (assuming 
regular burst occurrence) fell within data gaps. Thus, it 
was still possible for the bursts to be recurring on a regular 
timescale. 

1 The response used for extrapolation covered the range 0.1-200 keV with 200 logarithmically spaced energy bins. 



In order to independently test for the presence of these 
intermediate bursts, we examined ASM dwells during 
these longer intervals. Very few dwells occurred close to 
the predicted time of the bursts. Thus, we found no ev- 
idence for bursts at the predicted intermediate times; on 
the other hand, with the present data we cannot exclude 
the presence of bursts within these gaps. Taking into ac- 
count these missed bursts, we conclude that the observed 
burst times are consistent with regular recurrence inter- 
vals of 5.74 ± 0.13 hr in 1997-98, 4.10 ± 0.08 hr in 2000, 
and 3.56 ± 0.03 hr in 2002. 

When combined with BeppoSAX and ASCA observa- 
tions between 1996-98, as well as more recent observa- 
tions with RXTE, these measurements provide evidence 
for a long-term, steady decrease in At with time (Fig. 
3). Over the same interval, the 10-d averaged 2-10 keV 
ASM intensity measurements of GS 1826—24 show evi- 
dence for a long-term trend of increasing persistent flux. 
This trend was confirmed by the PCA measurements. The 
2.5-25 keV persistent flux measured during 1997-98 was 
(1.32 ± 0.05) x 10~ 9 ergs cm -2 s -1 , while in 2000 it was 
(1.80 ± 0.05) x 10~ 9 ergscm _2 s _1 . This variation is sig- 
nificant at the 6.4cr level. The most recent observation, on 
2002 July 29, found the persistent 2.5-25 keV flux to be 
higher again, at (2.19 ± 0.02) x 10~ 9 ergscm~ 2 s _1 . 

3.2. Variation of At with F p 

The burst interval measured from the RXTE observa- 
tions decreased by 40% between 1997-98 and 2002, while 
over the same period the mean persistent flux increased 
by 66%. In order to estimate the bolometric correction to 
apply to the 2.5-25 keV flux, and thus estimate the ratio 
a between the integrated persistent flux to the burst flu- 
ence (equation 1), we extracted PCA and HEXTE spec- 
tra covering the entire observation from five representa- 
tive pointings (1997 November, 1998 June, 2000 July, 
2000 September and 2002 July). We fitted each spec- 
tra to a model consisting of a Comptonisation compo- 
nent ("compTT" in xspec) and a Gaussian, both atten- 
uated by neutral absorption along the line of sight. The 
compTT component features a cutoff at around 2.8 times 
the temperature of the scattering plasma, which we mea- 
sured as 13.4 ± 0.5 keV in the mean. In each case we 
obtained a good fit (mean x 2 /dof = 1.02 ± 0.13) over the 
energy range 2.5-60 keV. The mean neutral column den- 
sity was n H = (1.52 ± 0.25) x 10 22 cm~ 2 , which is con- 
sistent with the mean value derived from the PCA-only 
pre-burst persistent emission fits. We then measured the 
(absorbed) flux in the 2.5-25 keV band, and compared 
this to the flux in the 0.1-200 keV band estimated using 
an idealized response 1 in order to extrapolate outside the 
PCA/HEXTE bands. We found the bolometric correction 
to be 1.678±0.016, where the error represents the standard 
deviation from the results for the 5 observations. 

The resulting variation of the burst recurrence interval, 
as a function of the bolometric persistent flux, is shown 
in Fig. 4; qualitatively similar results were found by Cor- 
nelisse et al. (2003) from BeppoSAX WFC observations. 
For the errors on the individual At measurements from 
RXTE, we adopted a fractional error corresponding to the 
standard deviation measured by Cocchi et al. (2001) of 
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Fig. 4. — Variation of the burst recurrence time {upper panel) 
and the burst fluence (lower panel) as a function of the esti- 
mated bolometric persistent flux in GS 1826—24, from RXTE 
measurements between 1997-2002. Error bars indicate the la 
errors. The curves show theoretical calculations for a range of 
metallicities: Z = 0.02,0.01,0.003, and 0.001. The solid angle 
(R/d) and gravitational energy have been chosen in each case to 



match the observed fluence and recurrence time at F v 



2.25 X 



For Z = 0.02,0.01,0.003, and 0.001, this gives 
R/d = 13, 10, 8, 6 km @ 10 kpc, and Q grav = 175, 196, 211, 215 MeV 
per nucleon. 

0.1 hr, for a mean delay time of 5.7 hr. This agrees well 
with the scatter in At we observe within each epoch. For 
cases where we m iss n o r more intervening bursts, we scale 
the error by 1/y/n + 1, since we assume that we are ob- 
serving n + 1 burst intervals in total. From a least-squares 
fit to At as a function of F p , using At oc F£ we found a 
power law index of 7 = —1.05 ± 0.02. 

From the estimated bolometric fluxes we then calculated 
the ratio a = E p /E^ from each pair of bursts measured by 
RXTE (Fig. 5). We found that a varied significantly with 
epoch, decreasing from w 44 to w 40 between 1997-98 and 
2002, as the persistent flux increased. The weighted mean 
was 41.7 ± 1.6. 

From the BeppoSAX/WFC observations, Ubertini et al. 
(1999) derived a value of 60 ± 7. Roughly consistent val- 
ues were determined from BeppoSAX/NFl, RXTE, and 
ASCA observations, of between 50-54 (in 't Zand et al. 
1999; Kong et al. 2000). To understand the discrepancy 
between these measurements and the value of a calculated 
in the present study, we note firstly that the previous mea- 
surements of the burst fluence from ASCA and RXTE by 



Kong et al. (2000) agree well with our measurements. Fur- 
thermore, all the persistent flux measurements appear to 
agree in general once the difference in the energy band is 
taken into account. However, the burst fluence measured 
by both the BeppoSAX instruments are around 40% lower 
than for ASCA and RXTE. The RXTE/VCk is known to 
measure fluxes that are systematically ~ 20% higher than 
some other instruments (e.g. Kuulkcrs et al. 2003); how- 
ever, this offset is insufficient to explain the discrepancy 
in the measured fluence. Furthermore, the estimate of 
a from RXTE should be independent of any systematic 
flux offset. Despite the substantially higher burst fluence 
from the RXTE and ASCA measurements of Kong et al. 
(2000), the calculated a was still close to that measured by 
BeppoSAX (Ubertini et al. 1999; in 't Zand et al. 1999). 
This appears to result from the bolometric correction on 
the 2-10 keV RXTE/ ASCA flux, which is not quoted in 
the Kong et al. (2000) paper but we estimate at 4-6. By 
comparison, the bolometric correction implied by the 2— 
10 keV and 0.1-200 keV BeppoSAX/NFl measurements 
by in 't Zand et al. (1999) is 3.3. From our broadband 
spectral fits in §3.2, we estimate a bolometric correction 
for RXTE flux in the 2-10 keV band as 3.06 ±0.02. Thus, 
we attribute the higher a measured by Kong et al. (2000) 
to an excessive bolometric correction factor. 



4. COMPARISON TO THEORETICAL IGNITION MODELS 

In this section, we compare the observed burst prop- 
erties with theoretical models of type I burst ignition. 
We calculate ignition conditions following Cumming & 
Bildsten (2000), and refer the reader to that paper for 
details. Since the calculation depends only on the lo- 
cal vertical structure of the layer, we give the results in 
terms of the local accretion rate per unit area m, and 
the mass per unit area or column depth y. We assume 
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Fig. 5. — Ratio of persistent to burst luminosity a = L p /L b (equa- 
tion 1), calculated from RXTE observations between 1997 and 2002. 
Error bars represent the estimated la uncertainties. The curves 
show theoretical calculations for the same values of metallicity as in 
Fig. 4. 
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a 1.4 Mq neutron star with radius R = 10 km, giving a 
surface gravity g = {GM/R 2 )(l + z) = 2.45 x 10 14 cm s~ 2 , 
where 1 + z = (1 - 2GM / Re 2 )- 1 ' 2 = 1.31 is the gravita- 
tional rcdshift. This value for the redshift is close to that 
recently measured for EXO 0748-676 (z ~ 0.35; Cottam 
et al. 2002). 

We calculate the temperature profile of the accumulat- 
ing layer of hydrogen and helium, and adjust its thickness 
until a thermal runaway occurs at the base. The temper- 
ature is mostly set by hydrogen burning via the hot CNO 
cycle, and therefore the CNO mass fraction Z, which we 
refer to as the metallicity. Our models also include com- 
pressional heating and a flux from the crust F crust , but the 
results are not sensitive to these contributions. A factor of 
two change in -F C mst gives a 2% (25%) change in ignition 
depth and burst energy for Z = 0.02 (Z = 0.001), with 
a much smaller change in the trend of these properties 
with to. We take -Fcrust to be constant over the timescale 
of the observations, i.e. F crust = (to) Q C rust, where the 
time-averaged local accretion rate (to) is set equal to the 
value for which the burst recurrence time is 5.7 hr, and 
Qcrust =0.1 MeV per nucleon (Brown 2000). 

To calculate the burst energy, we assume complete burn- 
ing of the H/He fuel layer, and that the accreted material 
covers the whole surface of the star. The total energy is 
then 4:TrR 2 yQ nuc /(1+z), where y is the ignition column 
depth, Q nuc is the energy per gram from nuclear burning, 
and the factor £b accounts for anisotropic burst emission. 
We set the latter parameter equal to unity initially, but 
use it to track where burst anisotropy might play a role. 
We write Qnuc = 1.6 + 4 (X) MeV per nucleon, where (X) 
is the mass-weighted mean hydrogen fraction at ignition. 
This expression for Q nU c assumes w 35% energy loss due 
to neutrinos during the rp process (e.g. Fujimoto et al. 
1987), and gives 4.4 MeV per nucleon for approximately 
solar hydrogen abundance (X — 0.7). 

Table 1 gives ignition conditions for four different metal- 
licities. In each case, we adjust to until At = (y/m)(l + 
z) = 5.7 hr, as observed in 1997-98. Lower metallicity 
models have less hot CNO heating during accumulation, 
giving a larger ignition column and therefore larger to to 
match At = 5.7 hr. We are free to vary to because of 
the uncertainty in the F p -m relation. The ratio of ob- 
served flux and inferred to constrains a combination of 
the gravitational energy release per nucleon, and the de- 
gree of anisotropy in the persistent or burst emission. The 
persistent flux is 

where £ p accounts for anisotropy in the persistent emis- 
sion, and Qgrav = c 2 z/(l + z) w GM/R'is the gravitational 
energy release per gram. The burst fluence is 

Inserting equations (2) and (3) into a = F p At/Eb gives 

Q Q„„=(|)A = 2 9 MeV( 5 i I )(|). ,4) 

The last column in Table 1 gives aQ nuc for the mod- 
els. Reconciling these values with equation (4) requires 



Cb/Cp = 0.5-0.6, within the range of expected values 
(e.g. Lapidus & Sunyaev 1985; Fujimoto 1988), especially 
since the "anisotropy" factor for the persistent emission 
includes other effects such as energy radiated outside the 
X-ray band. 

The predicted variations in fluence, recurrence time, and 
a with M are compared to the observations in Figures 4 
and 5. We adjust each curve to match the observed fluence, 
recurrence time, and a at F p — 2.25 x 10~ 9 erg cm~ 2 s _1 . 

The model with Z = 0.02 (approximately equivalent to 
solar metallicity) agrees very well with the observed burst 
energy and variation in a with M. The predicted burst 
energy (Table 1) is 5.3 x 10 39 ergs, corresponding almost 
exactly to the observed fluence for a distance of 6 kpc. The 
w 10% variation in a arises because hydrogen burning oc- 
curs between bursts as the layer accumulates. In a given 
fluid element, the time to burn all the hydrogen via the 
beta-limited hot CNO cycle is 11 hr (Z/0.02)- 1 (V /0.7) 
(e.g. Bildsten 1998), where Z is the CNO mass fraction, 
and X is the accreted hydrogen fraction. This means 
that as the recurrence time decreases from 5.7 to 4 hr, 
significantly more hydrogen per unit mass is present when 
the burst ignites, leading to an increase in <5 nuc , and cor- 
responding decrease in a. Although our models do not 
address the time evolution of the burst, the variation in 
mean burst lightcurves shown in Fig. 2 is consistent with 
a change in the fuel composition. 

However, the solar metallicity model is less successful at 
explaining the observed variations in At and E\, (Fig. 4). 
Whereas the observations show a decreasing ignition mass, 
At oc M -105 , and a w 3% rms variation in E^,, the models 
have At decreasing less steeply than 1/M, and a sa 10% 
increase in Ey, with M. These variations arise because 
as recurrence time drops, less hydrogen burning occurs, 
less helium is present at ignition, and the ignition mass 
increases. 

Extra heating of the accumulating layer would give a 
reduced ignition mass and better agreement with the ob- 
served trends of At and E\>. For the Z = 0.02 model in 
Table 1, the extra flux required to match the observed 
change in At is w 0.5 MeV per nucleon for an increase of 
50% in M. This could perhaps be provided by residual 
heat from the ashes of previous bursts (Taam et al. 1993; 
Woosley et al. 2003), although time-dependent simulations 
are required to test this. 

A second possibility is that the fraction of the neutron 
star surface covered by accreted fuel changes with M, as 
suggested by Bildsten (2000). We can infer the covering 
fraction from measurements of the blackbody normaliza- 
tion in the tail of the bursts. We show in Fig. 6 the mean 
value of i?bb/rfiok P c between w 20-50 s after the start of 
each burst. We staggered the time window over which each 
average was calculated, to ensure that the mean blackbody 
temperature over each time interval was approximately 
constant. We found a significant decrease of w 20% in 
the blackbody normalization. A variation of the covering 
fraction by this amount, and the resulting increase in lo- 
cal accretion rate to, would go some way to bringing the 
At and E\> variations into better agreement. However, the 
detailed relationship between the measured i?bb and the 
covering fraction is uncertain. 
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Fig. 6. — Blackbody normalization -Rbb/^lOkpc as a function of 
bolometric persistent flux, for bursts from GS 1826—24. The nor- 
malization was averaged over a time window fs 20-50 s following 
the start of the burst. Error bars indicate the la uncertainties. 

While the observed variation in a argues strongly 
against low metallicity models (which have a constant 
composition at ignition since little hot CNO burning oc- 
curs during accumulation), such models give much better 
agreement with the observed trends in At and E^. The 
constant composition leads to an ignition mass almost in- 
dependent of M, giving a constant burst fluence and At 
variations close to l/M. Low metallicity models overpre- 
dict the burst energy, since the reduced CNO heating leads 
to a larger ignition mass. Accounting for this discrepancy 
by increasing the distance to the source requires d w 10- 
12 kpc, above the current upper limit. However, there 
are other possible ways to bring the burst energy into 
agreement. Recent time-dependent simulations of bursts 
by Woosley et al. (2003) show that leftover hydrogen and 
CNO nuclei in the ashes of the previous burst lead to ex- 
tra heating, reducing the critical mass. Also, the fuel may 
cover only part of the neutron star surface. Table 1 lists 
the value of £ b ' (R/d) that gives agreement with the ob- 
served fluence (eq. [3]), where we now interpret R as a mea- 
sure of the area covered by accreted material rather than 
the neutron star radius. For the low metallicity models, 
we find ^^R w 5-6 km at 6 kpc, or covering fractions 
(assuming a neutron star radius of 10 km) of 30-40%. 

Finally, we address the sensitivity of our results to the 
fraction of accreted hydrogen. Recent evolutionary models 
involving intermediate mass binaries (Podsiadlowski et al. 
2002) propose that the secondary has undergone some 
main sequence hydrogen burning, leading to a relatively 
hydrogen-poor donor (X w 0.35-0.7). As an example, 
we show results for an accreted hydrogen mass fraction 
Xq = 0.5 in the lower half of Table 1. The ignition condi- 
tions are not very sensitive to the total amount of hydro- 
gen, since only a small amount of hydrogen burns during 
accumulation. However, the burst energy is significantly 
reduced, since Q nU c depends strongly on the amount of 



hydrogen versus helium. The variation in burst properties 
with flux is not a useful way to distinguish the amount of 
hydrogen, and is almost identical with Figures 4 and 5 for 
the X = 0.5 models. 

5. DISCUSSION 

The regularity of the bursting from GS 1826—24 makes 
it a unique source to compare burst theory and observa- 
tions. In this paper, we have presented RXTE/PCA ob- 
servations of 24 type I X-ray bursts from GS 1826-238. 
For the first time for this source, we have been able to 
study the properties of bursts in detail at different mass 
accretion rates. The bursts from a given epoch have iden- 
tical lightcurves (Fig. 1 and 2), and show long tails likely 
powered by rp-process hydrogen burning, as suggested by 
Bildstcn (2000). The bursts from a given epoch are con- 
sistent with a single recurrence time, and this varies with 
accretion rate as At cx M -1 ' 05 . The burst fluence in- 
creased slightly with the accretion rate, and a decreased 
by w 10%, with mean value 41.7 ± 1.6. 

The decrease in a with M implies that stable burning 
of hydrogen occurs between bursts. This requires solar 
metallicity (Z « 0.02) in the accreted layer (Fig. 5). Solar 
metallicity models also give good agreement with the ob- 
served burst energies, as pointed out by Bildsten (2000). 
However, the small variation in fluence, and At variation 
with M is difficult to explain with solar metallicity mod- 
els, which predict 10% variations in fluence, and At de- 
creasing less steeply than l/M. Low metallicity models 
(Z ~ 0.001) naturally explain the small variation in burst 
fluence, and w l/M scaling of At, since little CNO burn- 
ing occurs between bursts giving an ignition mass almost 
independent of M. However, the observed a variations 
with M rule out low metallicity models. 

There are several possible ways to reconcile the solar 
metallicity models with the observations of At and E\>. 
Studies of timing and spectra of LMXBs indicate that 
Lx is not always a good indicator of M (e.g. van der 
Klis et al. 1990), we have assumed Lx oc M here. Ex- 
tra heating of the accumulating layer would act to re- 
duce the critical mass and bring the observations and the- 
ory into agreement. One possibility is that residual heat 
from the ashes of previous bursts heats the layer (Taam et 
al. 1993; Woosley et al. 2003), time-dependent simulations 
are needed to test this. If the fraction of the neutron star 
surface covered by fuel changes with M, the changing local 
accretion rate per unit area could also reconcile the mod- 
els and observations. We find that the blackbody radius 
i?bb in the tail of the bursts decreased by «20% between 
the observed epochs. If this indicates a change in covering 
fraction, it would almost be enough to explain the discrep- 
ancy. However, the covering fraction decrease with M is 
opposite to the increase suggested by Bildstcn (2000) to 
explain trends in burst properties. Furthermore, i?bb is at 
best an uncertain measure of the emitting area. Due to 
variations in the atmospheric opacity with photon energy, 
the measured blackbody temperature is generally higher 
than the effective temperature (e.g. London et al. 1986; 
Shaposhnikov et al. 2003). Consequently, i?bb is expected 
to give an underestimate of the emitting area, by a factor 
which depends on the detailed structure and composition 
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of the NS atmosphere. Better R/d measurements could 
be obtained with high resolution spectra and appropriate 
model atmospheres. 

Woosley et al. (2003) have recently calculated time- 
dependent multi-zone burst models including a large nu- 
clear reaction network. They show that both the rise 
and decay times of the burst depend sensitively on nu- 
clear properties in the rp-process path, and in addition 
find burst rise times of several seconds, in good agreement 
with the observations. There is much to learn from a de- 
tailed comparison of time-dependent simulations with the 
burst lightcurves. In particular, the variation in decay 
time with accretion rate may directly reflect the changing 
abundance of hydrogen at ignition. 

Many bursters show a transition from regular, frequent 
bursts with long durations at low to, to irregular, infre- 
quent and short bursts at high to, a behavior first dis- 
covered by EXOSAT (van Paradijs et al. 1988). By com- 
paring BeppoSAX/MVFC observations of nine regular type 
I bursters, Cornelisse et al. (2003) propose that this be- 
havior is common to all bursters, and that the transition 
occurs at a universal luminosity Lx « 2 x 10 37 erg s _1 . 
While the physics of the transition to irregular bursting 
is not understood (e.g. see discussion by Bildsten 2000), 
we have shown in this paper that the regular bursting 
seen in GS 1826—24 is well-understood as being due to 



helium ignition in a hydrogen- rich environment. This ar- 
gues that other bursters such as KS 1731-260 are in this 
burning regime when regular bursting is seen (and not 
bursting via unstable hydrogen ignition, as suggested by 
van Paradijs et al. 1988 and Cornelisse et al. 2003). Bep- 
poSAX found a peak burst rate for KS 1731-260 (the 
burster with the largest range in observed to's) of w 9 
bursts per day. Extrapolating the observed 1/M de- 
pendence of At, we expect GS 1826—24 to reach this 
burst rate, and perhaps transition to irregular bursting 
at F p « 4.6 x 10~ 9 erg cm~ 2 s _1 . Further monitoring of 
the burst behavior of GS 1826—24 will hopefully probe the 
physics of the transition to irregular bursting. 
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Table 1 

Ignition conditions for bursts with At = 5.7 hours' 



z 


ni4 


ys 


Ts 


X (X) Qnuc b 




R/d d 














(10 39 erg) 
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Accreted H fraction 


X„ = 0.7 






0.02 


0.97 


1.5 


2.2 


0.40 0.55 3.8 


5.3 


9.6 


152 


0.01 


1.4 


2.1 


2.1 


0.55 0.62 4.1 


8.0 


7.8 


164 


0.003 


2.1 


3.3 


2.0 


0.65 0.68 4.3 


13 


6.1 


172 


0.001 


2.8 


4.3 


1.9 


0.68 0.69 4.4 


17 


5.3 


175 










Accreted H traction 


X„ = 0.5 






0.02 


0.92 


1.4 


2.1 


0.20 0.35 3.0 


4.0 


11 


120 


0.01 


1.2 


1.9 


2.0 


0.35 0.42 3.3 


5.8 


9.2 


132 


0.003 


1.8 


2.9 


1.9 


0.45 0.48 3.5 


9.3 


7.3 


140 


0.001 


2.4 


3.8 


1.8 


0.48 0.49 3.6 


12 


6.3 


143 



Note. — The first six columns arc: CNO mass fraction Z\ local rest mass accretion rate at 
the surface of the star rh.4 — m/10 4 g cm -2 s _1 ; ignition column depth y§ — y/10 8 g cm -2 ; 
ignition temperature T s — T/10 8 K; hydrogen mass fraction at the base X; mean hydrogen mass 
fraction (X) — f dyX(y) fy. 

a Ai — (y/m)(l + z) is the recurrence time measured by an observer at infinity. 

b Nuclcar energy release, Q nuc — 1-6 + 4.0 (X) MeV per nuclcon. 

c Predictcd burst energy assuming the fuel covers the surface of the star, -EWrst — 
/(l + z) (we take R = 10 km, z — 0.31, £b — 1)- The observed burst energy 
is 4.9 x 10 39 erg (d/6 kpc) 2 . 

d The R/d that gives agreement with the observed burst fluence. This quantity is £~ 1,/2 (Rfd) 
if the emission is anisotropic. 

e We adopt the observed value a — 40. 



